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Contributions to the Study of Oscillation Properties of the 

Solutions of Linear Differential Equations 

of the Second Order.* 

By R. Gr. D. Eichaedson. 



Introduction. 

The study of boundary problems for linear differential equations of the 
second order dates back to the time of Euler and D'Alembert, with whom it 
arose in connection with problems of mathematical physics. Beginning with 
the fundamental paper of Sturm in 1836, there have been extensive investiga- 
tions t in this field in recent years, notably by Klein, Bocher, Stekeloff, Kneser, 
Hilbert and Birkhoff. Since the differential equation of the second order is 
of such fundamental importance in so many fields, and since similar general 
problems for equations of higher order can not be handled by processes so far 
devised, the invention of new methods and further investigation of the nature 
of solutions find ready justification. 

The equation to be studied will be taken in the form 

^(p{*)t^) + G(*,Vv(*)-(W.).+G(*,*)V=0, 0<x<l, (1) 

where G{x, X) is a function depending on a parameter %. The solution y(x) 
of this self-adjoint equation shall be subject to the self-adjoint boundary con- 

diti ° nS «tf (0) +«#.«)) +032/(1) +«#„(!) =0, ' 

fty(0)+fty„(0) +&</(!) +&y„(i)=o, (2) 

p(l) (ajft-aA) =p(0) (a 3 /3 4 — a 4 /? 3 ), . 

where the two sets of real coefficients a 2 , a 2 , a 3 , a 4 and ft , /? 2 , ft > ft are 
linearly independent. The most important special cases of these boundary 
conditions are given by 
y(0)=y(l)=0; y(0)=y.(l)=0; y.(0) =y(l) =0; y x (0)=y x (l)=0. (2') 

* Read before the American Mathematical Society, September 4, 1917. 

t Tor existing methods and literature of the subject see B6cher, EncyhlopSdie Mathematischen 
Wissenschaften, II A7a, Proceedings International Congress of Mathematicians, Vol. I (1912), p. 163, and 
" Lecsns sur les Methodes de Sturm " (1917) ; Lichtenstein, Rendiconti del Circolo Matetnatico di Palermo, 
Vol. XXXVIII, p. 113. 
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There are many interesting questions in regard to this linear problem. Do 
there exist parameter values a such that there are solutions of (1) satisfying 
relations (2)? If so, how many are there, and how are they distributed? What 
is the nature of the corresponding solutions? How do the solutions vary with 
change of the coefficients of the equation and of the boundary conditions? 
When does the totality of solutions form a fundamental set in terms of which 
functions may be expanded? 

Among the methods which have been used in studying the problem are : 

(1) Differential equations including the use of comparison, approximation and 
asymptotic expressions ; (2) the minimum principle in the calculus of varia- 
tions; (3) integral equations; (4) the theory of linear algebraic equations 
in an infinite number of variables;* (5) a limiting process with linear 
algebraic approximating difference equations (cf. § 1, IV). The methods 
developed in this memoir would fall under (1) and (2) and may be character- 
ized as a free use of the differentiation of fundamental formulae with regard 
to the parameters involved in the equations and in the boundary conditions, 
together with the proof that under the conditions imposed certain integrals 
are positive. 

Exact oscillation theorems for solutions under the boundary conditions 

(2) have been developed by Birkhoff f for the case of the special equation 

7)C 
y xx +G(x,7i)y = 0, ~>0, lim <?=-», lim <?= + «>. 

OA \=— co X=+oo 

Another important special case of equation (1) in which 6(x, a) contains the 
parameter linearly 

(P0,).+ (2+M)y=O, (3) 

has been studied very extensively. When k (x) > this equation may be 
reduced to a form included in that investigated by Birkhoff. 
The definite case of (3), viz., when one of the integrals 



C ky*dx, f (pyl—qy 2 )dx 



has one sign for all functions y(x) considered, has been discussed in many 
phases by mathematicians since the time of Sturm. By means of his theory 
of integral equations Hilbert $ established the existence of characteristic 

* Lichtenstein, loc. cit. 

f Transactions of the American Mathematical Society, Vol. X (1909), p. 259. It should be remarked 
that both in the second and third lines from the end of the statement of the principal theorem (p. 269) 
instead of p-f- 1 we should read p — 1. 

$ " Grundziige einer allgemeinen Theorie der linearen Integralgleichungen " (Teubner, 1912) . 
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parameter values and characteristic solutions of equation (3) for this definite 
case. As will be shown in § 5 the boundary conditions which he used are nor- 
mal forms of (2). Oscillation theorems for equation (3), under some of the 
simple boundary conditions (2'), have been established by various means; 
among others by setting up the corresponding calculus of variations problem, 
and interpreting the Jacobi criterion.* 

When k has both signs and q is positive in at least a part of the interval 
and sufficiently large there, it is not necessary that either of the integrals in 
question be definite. This case, which we shall call the non-definite, was first 
discussed incidentally by the author in a paper f in which he was treating the 
problem of oscillation theorems for two equations with two parameters. He 
showed that when the boundary conditions are «/(0) =t/(l) =0, there exists an 
integer % such that for n<n x there are no real solutions which have n zeros, 
while for n>n 1 there are at least two.\ At that time all the principal results 
of §§ 2-4 were obtained, but were not published. 

For the special equation (3) and k>0 the theorems proved by Birkhoff 
were rediscovered by Haupt§ in his dissertation. The oscillation theorem 
stated in this dissertation for the case that h changes sign is corrected in a 
later article,|| and by using the methods which I had developed, various oscil- 
lation theorems for the general equation (1) are derived. It is also shown by 
means of expansion theorems that if the non-definite equation (3) be taken in 
a certain normal form there exists an integer n 2 such that for n>n 2 there are 
precisely two solutions with n zeros and satisfying the boundary conditions (2). 

The object of the present memoir is to investigate the conditions to be 
imposed on G(x, X) regarded as a function of %, so that definite oscillation 
theorems for solution of (1) may be determined. In § 2 an attempt is made 
to bring as close together as possible necessary conditions and sufficient condi- 
tions for a limited or unlimited number of oscillations. A criterion for the 
behavior of the zeros with change of the parameter is obtained in §3. This 
permits the development of very general theorems for the unique existence of 
solutions vanishing at the end points and possessing a prescribed number of 
zeros, and also for the existence of two, and only two solutions of this nature. 

* Mathematisehe Annalen, Vol. LXVIII (1910), p. 279. 

f Transactions of the American Mathematical Society, Vol. XIII (1912) , p. 22. 

% That there were exactly two when »>n 1 was stated by the author in a paper in the Mathematisehe 
Annalen, Vol. LXXIII, p. 289, in which he was discussing oscillation theorems for three linear equations 
with three parameters. The subsequent results of the memoir were not affected by this error which was 
corrected by a note in Vol. LXXIV (1913), p. 312, of the same journal. 

§ " Untersuchungen ttber Oszillationstheoreme " (Teubner, 1911). 

|| Haupt, Mathematisehe Annalen, Vol. LXXV1, p. 67. 
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These theorems contain as special cases all known results in this field, and 
some new special cases are set forth in detail. 

The non-definite case of (3) and (1) is discussed in §4. The question of 
whether there may be for a given oscillation number more than two parameter 
values for which there are solutions of (3) (2) is settled by giving an example 
in which for any n in an interval n lf « 2 there are four values of Jl corre- 
sponding. That there exists an integer n 2 such that for n>n 2 , there are 
exactly two solutions is proved by a method entirely different from that of 
Haupt, and in some particulars it would seem that the resulting theorem is 
less satisfactory, in others more satisfactory than his. The theory is also 
extended to cover some corresponding cases of (1). Concerning the complex 
solutions which correspond to values of n < n x some theorems are derived. 

In the latter half of the memoir a method is developed for obtaining the 
facts in regard to the solutions of the general equation (1) under the general 
boundary conditions (2). With this end in view § 5 is devoted to a reduction 
of the boundary conditions to normal forms by means of the usual transfor- 
mation of the dependent variable y, which leaves the number of zeros 
unchanged. Bach of these three normal forms, 

I. ffy(0)+y.(0)=0, Ty(l)+y.(l)=0; 

II. y(0)=hy(l), hp(0)y,(0)=p{l)yM; 

III. y(0)=fo(l)y.(0), lp(0)y x (0)=y(l) 

contains one or two parameters a, t, h, I. When special values or oo are 
assigned to these parameters, the forms reduce to the simple cases (2') for 
which the facts are readily obtainable from the developments of the earlier 
sections. The first form is of essentially different character from the others. 
In the latter the parameter h or I is a double-valued function of \ which 
in general is real in sub-intervals only, while in the former each of the 
parameters c, t is a single-valued function of the other and of X, and is real 
throughout. By letting X vary, and calculating the rates of change of these 
parameters and of G(x, X), theorems of oscillation are derived (§§6-8) for 
boundary conditions in each of the normal forms. While detailed results are 
not given in all cases, these are immediate developments of the fundamental 
facts ascertained. 

Any linear differential equation of the second order 

*y«.+4tf.+fy=0 or W„+<te,+ (M 1 +d t )y=0, i»0, 

may be thrown into the corresponding self-adjoint form (1) or (3) on multi- 

plying by the function p(x)=eJ° * * ; the corresponding self -adjoint boundary 
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conditions (2) are changed in form only by the substitution p(0) — l, 

p(l)=e J <" 1 ' . 

The oscillation theorems remain unchanged if the variables are subjected 
to the usual transformations of the dependent and independent variables 

y=vjy, x=%(x) (y;(x)=f=0, -pnfcO ). The resulting equation is of the form 

fe+^+%=0 or ^y»+yg-,+ (^i+9 g )y=0, 
the new boundary conditions for the new interval x , x x being of the form 
ai2/(«o) +O2^(«o) +as^(^i) +a^ i (x 1 ) =0, 
Ay&o) +te(«o) +&#(«i) +fe(«i) =0, 

as may be shown by computation.* The invariantive property of self-adjoint- 
ness gives to the results obtained for (1) or (3) a very general character. 

§ 1. Some Properties of Solutions of the Differential Equation. 
In his original memoir Sturm studied the differential equation in the form 

where K>0 and G both depend on a parameter PI. But by a change of 
variables this may be reduced to the form 

(py.),+O(x,X)y = 0, (4) 

where p(x) is a positive function independent of /t. We lose nothing in 
generality by considering this latter equation. Regarded as functions of x 
the coefficients p, G will be postulated as continuous together with as many 
derivatives as is desired, while G will be considered as analytic with respect 
to /U The usual modifications of the results derived can be written down im- 
mediately if less stringent hypotheses are imposed. The trivial solution t/=0 
will be excluded from the discussion. Some theorems concerning solutions of 
(4) will now be reviewed. 

I. If one zero of a solution of (4) is held fixed, all others are moved 
nearer to it by a decrease of p or an increase of G. If, for example, p(x) is 
less than a constant P, and G(x, Pi) is greater than a constant y >0, the zeros 
of (4) are closer together than those of the equation 



* For the special case of a transformation of the dependent variable only, this computation is 
given in § 5. 
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which has a solution « = sin Vg(x4-c) with zeros at intervals of —7= . By 

fixing p, and taking G large enough in any interval of x, the zeros of (4) may 
then be made as close together as desired. 

II. The special equation where G{x,X) contains the parameter linearly, 

(py.).+ (g+a.ft)y=o, (5) 

has been much discussed. The boundary conditions 

y(0)=y(l)=0 (6) 

are of the greatest interest and three cases may be distinguished. 

(A) Orthogonal Case, -when k(x)>0. There is then an infinite number 
of parameter values ^ m {^i<^< • • • ■ •) with a limiting point at positive infinity 
only, for each of which a solution Y m satisfying (6) exists. The number of 
zeros of the solution Y m (including those at x=0 and a?=l) is m + 1. 

(B) Polar Case, when k(x) takes on both signs and the integral 

D{y)=C(pyl~<itf)dx (7) 

is positive-definite,* that is, for the given boundary condition (6) B(y) can 
not take on negative values. There are two sets each of an infinite number of 
parameter values < "k x < A, 2 . . . . , > 7*_ x > A_ 2 . . . . , with limiting points at 
positive and negative infinity respectively, corresponding to which solutions 
Y m , r_ m exist. For both Y m and Y_ m the number of zeros is m+1. 

(C) Non-definite Case, when both of the integrals 

f ky 2 dx, f (pyl—qy 2 )dx 

may take on negative values. This will be discussed in § 4. 

III. For the equation (5) and the special boundary conditions (6) cer- 
tain minimum properties may be stated. In the orthogonal case the minimum 
of the integral D(y) for those values of y which satisfy (6), and the normal- 
izing and orthogonalizing conditions 

f ky 2 dx=l, I kYiydx=0, i — l,2,....,m — 1, (8) 

is l m , and is furnished by the normalized solution Y m of (5). In the polar 

case the minimum of D(y) for those values which satisfy conditions (6) and 

(8) is % m and is furnished by Y m ; the minimum subject to the conditions (6) 

and , j 

j ky 2 dx — — 1, ( kY_iydx—0, i — \, 2, . . . ., m — 1, 

* Bocher has pointed out (Proceedings International Congress of Mathematics, loc. cit., p. 173) that 
the special case of the polar problem where q <, can be treated by the method of Sturm. This remark, 
however, does not apply to the polar case in its most general form. 
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is — %_ m and is furnished by Y_ m . In the non-definite case, D (y) can be nega- 
tive and the minimum (even in the simplest problem (w = l)) may not exist. 
However, in view of the developments of § 4, it would seem probable that by a 
modification of the discussion, the solution Y m , for m large enough, may be 
regarded as furnishing a minimum of a calculus of variations problem. 

IV. While the notion of regarding a differential equation directly as the 
limit of a set of difference equations has been used heuristically since the time 
of Euler, it was, so far as the author is aware, first made definite and rigorous 
in a recent paper.* In that paper the problem actually discussed is that of 
existence theorems for partial diff erenl ial equations with given boundary con- 
ditions. But the same method applies, for example, to equation (1), and the 
discussion is essentially simpler for the equation in one dimension than for 
that in two or more. For the sake of simplicity let us confine ourselves to the 
case (5) (6) and consider the unit interval to be divided into m equal parts, 

the values of y, p, q, k at the point — (*=0, 1, , m) to be denoted by y t , p { , 

q i} k i} and difference equations 

m 2 [p i+1 (y i+1 — y { ) —Piiyi—Vi-i) ] + q i y i +M i y i =0 (»=1, 2, , m— 1), (9) 

to be set up. In order that there be solutions of these equations "k must be 
one of the m — 1 roots of the determinant formed from the coefficients. With 
increase of m the number of points at which y is defined increases, but we can 
pick out corresponding parameter values and solutions of the various sets of 
difference equations, and if proper continuity conditions are imposed on the 
coefficients of the differential equations, it may be shown that the correspond- 
ing sets of parameter values approach as a limit a parameter value of (5), and 
corresponding solutions approach a solution of (5). In this way the infinite 
set of solutions of the differential equation is obtained. If k has both signs 
and q is positive and sufficiently large, at least in some part of the interval, 
some of the Vs and the corresponding y's will be complex.! But in all those 
cases of equation (5) heretofore treated (the orthogonal and polar), the 
method of passing to the limit in (9) suffices and gives a simple proof of the 
fundamental facts. This method can be extended to a treatment of existence 
theorems for solutions of the equation (4) with more general boundary con- 
ditions. 

* Transactions of the American Mathematical Society, Vol. XVIII (1917), p. 489. 
f One method of proving this would be by noting that the parameter values and solutions are 
approximations to those of the differential equation which are shown in §4 to be complex. 

37 
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§ 2. Sufficient Conditions for the Existence of Solutions with an Unlimited 

Number of Zeros. 

In the consideration of the equation 

(py.).+G(x,Vy=o (10) 

let us impose the restriction that for finite values of X the function Q (x, X) is 
limited, and in particular that G(cc, 0) is limited. The interesting cases will 
be covered by one of two hypotheses, which will be justified by the later devel- 
opments of this section. 

Hypothesis A. For at least one point of the interval, the upper limit of 
G{x,"k) becomes infinite with X (Urn G(x, X) = + oo), and in such a manner 

X=+oo 

that the number of zeros of the solutions increases without limit with /I. 

dC 
The problem treated by Birkhoff* where -~r- >0 and G(x, +oo) = + oo is 

a special case, and the orthogonal problem (§1, II) is still more special. 

Hypothesis B. For at least one point of the interval Urn G(x,7i) = -\-oo; 

\=+oo 

for at least one other, Urn G(x, /I) = -f oo; and in both cases G(x, Jl) increases 

X=— 00 

in such a manner that the number of zeros of the solutions increases without 
limit with X. 

The polar case (§1, II) is included in this hypothesis. 

Theoeem I. In order that there be a set of parameter values X such that 
the number of oscillations of the corresponding solutions of (10) be unlimited, 
it is necessary that in the neighborhood of at least one point Urn G(x, h) = + oo 

X= + oo 

or Urn G (x, "k) = + °°. 

X = — 00 

For, as the number of zeros increases, the length of the smallest interval 
decreases without limit. The zeros of the equation 

Py xx -\-^{ x i %)y=0, P=maximumf>(a;) (11) 

are farther apart than those of (10) (§1, I). To establish the theorem for 
(10) it is then only necessary to prove it for (11). Let us denote by a, /? that 
pair of consecutive zeros of (11) whose distance is a minimum. In such an 
interval a, /3 [y(a) =y(8) =0], y may be taken positive, and since the equa- 
tion is homogeneous, its solution for all values of /I may be multiplied by a 
constant so that the maximum is 1. The maximum of y x must be at least as 

great as -^ which is the slope of the line joining (a, 0), (/?, 1). To investi- 

* Loo. dt. 
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gate the maximum of — — we note that y x is zero, at some point of the inter- 
val and at least as great as -~ at another. Hence — y xx must be at least as 

great as y-5 ^ at some point, and further max f — — )> ~_ 2 . Since 

for at least one sub-interval the length approaches zero, and since 

y 

it follows that for the interval 0, 1 lim max G (x, X) = + oo. Since by hypothesis 
G can become infinite only for X=±oo the theorem may be readily deduced 
from these results by the usual processes of reasoning. 

To show that the necessary condition of Theorem I is not sufficient, let us 
consider the following example : In the interval 0, £ — s we set up the function 

y= — sin -rr- , and in the interval i + e, 1 the function y — — sin -^ (1 — x) ; 

71 2 " ' 7t Z 

both of these arcs satisfy the equation y xx + -j- y = 0. If continued, they 

would meet at the point x=i with an angle arc tan 2\/2. If in the interval 
i — s, 1 + e an analytic curve is introduced which is tangent to these two arcs 

(on one of which y x > —^ and on the other y x < ~\, it must have at some 

point — y xx > ■ ■, and if e is taken small enough, max ( — — J is great at 

pleasure. Hence denoting by y xx +G(x, %)y=0 the differential equation which 
has for solution the function defined for the interval 0, 1 by this method, and 

setting %= — , we know that lim max G(x, /I) in the neighborhood of oc=i 

e e=0 

increases without limit. On the other hand none of the suite of functions y 
has any zeros in the interval. 

Having shown by this example that the necessary condition of Theorem I is 
not sufficient, let us now deduce a sufficient criterion that when lim G(x,X) = ca, 

the number of zeros of solutions of (10) be unlimited. A similar discussion 
can be given for the case lim G{x, X) — 00. Denoting as before by P the maxi- 

mum of p (x) , and by M an arbitrarily large constant, it is possible to find a X 
and an ^-interval such that G(x,X) >MP in this interval, whose length will be 
denoted by e M . The number of zeros of the solution sin -\/M(x+c) of the 
equation y ,,,.+ My =0 in an interval of length e M is not less than the integral 
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part of — — - . It follows immediately from § 1, I that the number of zeros of 
it 

solutions of (10) is not less than that of sin -\JM(x-\-c). Hence the 

Theorem II. In order that there be N zeros in a solution of (10) it is 
sufficient that one can find an M such that — — - >N, where e M is the length of 
a sub-interval where G(x, X) >MP. 

Corollary I. In order that the number of zeros be unlimited it is 
sufficient that Urn e M '\/M=<x>. 

Corollary II. If throughout any sub-interval of fixed length the value 
of G{x,X) increases without limit, the number of zeros increases indefinitely. 

However, the sum of the sub-intervals in which G>M may remain above 
a positive constant without compelling the number of zeros to increase with M . 
For example, one can set up a function which has no zeros except at x=0 and 
x—1 and which, except in the neighborhood of these points, oscillates between 
y=i and y=Y, the number of oscillations increasing indefinitely with X. More- 
over, by taking the oscillations frequent enough one can have G(x, X) >MP in 
portions which total at least one-quarter (or any other proper fraction) of the 
interval, the great curvature downward in these sub-intervals being counter- 
balanced by curvature upward in those where G(x, X) has the opposite sign. 

On the other hand, we have seen that the number of zeros of the solution 
sin ^M(x-\-c) of the equation y xx +My=0 is not less than the integral part of 

— — - where e M is the length of the interval. If this interval is divided into 

71 

two parts £ 1} e 2 for which there are solutions sin -\/M(x-\-c x ), sin ^M(x+c 2 ) 
respectively, the number of zeros can not be reduced by more than one. For, 

the sum of the integral parts of — — - and — — £• cannot differ by more than 

n 7i 

one from the integral part of — — — — . It follows in the same way that if 

71 

the interval is divided into n x parts the number of zeros of the solutions can 
not be decreased by more than %-— 1. But in the intervals where G(x, X) >MP 
the zeros of the solutions of (10) must be at least as many as the minimum 
number for y xx +My = 0. 

Theorem III. In order that the solution of (10) have N zeros it is 

\I~M.Y! 

sufficient that one can find an M such that — — — — n 1 +l>N, where y\ m denotes 

7t 

G 

the sum of the lengths of the n x intervals in which 7— r- >ilf. 

maxp(x) 
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§3. Behavior of the Zeros with Monotone Change of 2,. 
Let us fix the zero at the left-hand end of the interval 0, 1, denote such a 
solution by Y [Y(0) =0], and investigate what happens to the other zeros as 
a, increases. It is convenient to think of G(x, X.) as being defined for values 
x>l. Since the coefficients of (10) are continuous, the zeros move continu- 
ously, and since we can not have both y(x)=0 and y'(x)=0 without having 
t/=0, the zeros can not coalesce and then disappear. By differentiating (10) 
with regard to 2> we have 

and on multiplication of this equation by — y and of (10) by ~|, addition and 
integration from x x to x 2 we get the fundamental formula * 

[-i];-[-t]:=ri^- 

For the particular solution Y(x) we have J(0) =0, — ^j— -, and if a is another 
zero of Y the formula becomes 

p{a )YA«)^(«)=f o a ^Y>da. (13') 

If the integral on the right is positive, the signs of Y x (a) and — ^— - are the 

same. Hence when Y x (a) is negative, 7(a) is decreasing with increase of X, 
and this zero of Y is moving to the left; when Y x {a.) is positive, Y(a) is 
increasing and the zero is again moving to the left. In the same way it may 
it may be shown that all zeros for which the integral on the right of (13') is 
negative, move to the right. When the integral is zero, further investigation 
is needed. We have now proved 

*If p is also a function of X this formula, becomes 

In terms of this hypothesis results analogous to all the succeeding theorems of this paper may be written 
down. Since, however, by a transformation of variables p can be made independent of X we shall, for 
the sake of simplicity, confine ourselves to the problem proposed. 

Qv (0) Qy(0) 
If at <b = the condition y x (0) -f- ay (0) =0 is imposed where a is a constant, then — xf +<r Z. =0, 

and this formula (13) can be written 



■ yp to 



/®M\ r Qy Qy*l r % 90 2J 



9« 
Hence, in order that the roots of a solution y and the roots of ^ alternate it is sufficient that 



/. 



* 9<? 
-^r—y'dx have one sign. 

OA 
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Theorem IV. If Y(0) =0 and if a is a parameter value for which there 
are subsequent zeros of Y at a x ,...., a„ then with increase of a a zero, a { , 
moves to the left or right according as 

• a * dG{x, a) 



X 



is positive or negative. 

The argument used in deducing Theorem IV is still valid if in place of a 
zero of y we introduce a zero of y x . 

Theorem IV A. If y(0)=0 ory x (0)=Q and if a is a parameter value 

such that y(a) =0 or y x (a) —0, then with increase of a the zero a moves to the 

left or right according as 

dG(x,A) . 



S 



8a tf (*>*•)** 

is positive or negative. 

Corollary. In order that a be a multiple characteristic number and 
Y(x, a) a multiple solution of the problem (10) (6) it is necessary* that 



X 



ldG ^>T' ix ,X)da:=0. 



X 



G can be so chosen as a function of a that with increasing a the integral 

i dG 
, -5— Y z dx is alternately positive and negative, and since the corresponding 

OA 

intervals of a may be small at pleasure, we see that a zero may pass back and 
forth through x=l an infinite number of times. There may then be an infinite 
number of a's for which exist solutions of (10) (6) with a fixed number of 
zeros. 

As a runs through its interval there will be at least one value for which 
the number of zeros of Y{x, a) is a minimum. For both the orthogonal and 
polar cases this minimum number of oscillations is zero, but in general this 
will not be the case. If for a—W there is a minimum number n x of oscilla- 
tions, then from the principle of continuity we can argue that for n>n ± there 
is under Hypothesis A (§2) at least one value of a for which exists a solution 
of (10) vanishing at x=0 and x=l, and oscillating n times (having n zeros) 
in the interval, while under Hypothesis B there are at least two. 

One fundamental problem for investigation is the determination of 
sufficient conditions that in a given interval of a (which may include infinity) 
there is only one solution of (10) (6) oscillating a given number of times; in 

* This condition may be shown to be sufficient also. 
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other words, sufficient conditions that the zeros of Y(x) pass through the point 
x=l in one direction only as a increases (or decreases). Another problem is 
that of determining when there exist precisely two solutions oscillating a given 
number of times. 

Theorem V. // for values of a in an interval {which may include infinity) 
there are solutions of (10) (6) for which the minimum and maximum number 
of oscillations are n x and n 2 (n 2 may be infinite), respectively, and if the value 

dC 

of -*— is equal to or greater than $(a)G(x, a), where $ is a positive function 

of a, then there is one and only one solution oscillating n times {ni<n<n 2 ) in 
the interval 0<a?<l. 

For, we have on multiplying (10) by Y, and integrating under the boun- 
dary conditions (6), 

CpY\dx= CGYHx. (14) 

Jo Jo 

From the hypothesis we may then write 

J- 1 ?lG r 1 r 1 

^-YHx><p(A) f GY 2 dx=Q(A) f pY 2 x dx>0 
o OA Jo Jo 

and since this holds for all a, Theorem IV gives the desired result. 

The theorem just proved fits the case where with unlimited increase of a 
the maximum of G increases without limit. For the case where a approaches 
— qo, we observe from Theorem IV that a sufficient condition for a similar 

theorem is that J -^- Y 2 dx be negative. To insure this, it is sufficient to take 

dC 

■^r- < — ^G, where ^(a) is a positive function. Then 

J^Y 2 dx<—^(A)fGY 2 dx=-^(A)fpYldx<0. 

Combining this result with Theorem V we may state the following: 

Theorem VI. If for two values a', a" (a' > a") of A there are n x zeros of 

dC 

the solution Y within the interval 0, 1, and if for a>a', -^->$i(h)G(x, a), 

dC 
and for a,</l", -~— < —<p 2 (A)G(z, a) where $ 1} 4» 2 are positive functions defined 

in the intervals a', +oo and -co, a", respectively, then under Hypothesis B 
(§2) for »>»! there are two and only two values of a for which there is a 
solution Y(x,a) of (10) (6) with n zeros. 
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In the polar problem we have the special case of the hypothesis of this 

theorem where W=1" = 0, %=0, G=q + M, -g— =k, &= — >0, — &> = -=- <0. 

A more general special case of Theorem VI is discussed in § 4. 
Since formula (14) is valid also when 

y(0)=y x (l)=0, or y.(0)=y(l)=0, or y„(0)=y.(l)=0, 
theorems analogous to VI can be written down for any of these boundary 
conditions. 

§4. The Non-Definite Case. 
With regard to applications, the differential equation 

Ow.).+ (2+aA)y=o, (15) 

with boundary conditions 

2/(0) =«/(l) =0, (16) 

is the most important special type of (10). The problem where one of the 
integrals 

D(y, 0, 1) = f (pyl—qy 2 )dx, f V^ 

is definite, has been studied in great detail. Let us first investigate some 
properties of the remaining non-definite case where both these integrals may 
take on negative values and later deduce some analogous results for the more 
general equation (10). In §§6-8 these results will be extended to cover the 
more general boundary condition (2). 

For q<0 the integral D(y, 0, 1) can not be negative. If q is positive in 
at least a part of the interval 0, 1, the minimum of the integral $\py\d% under 
the conditions (16) and j*qy i di>=l is (§ 1, III) the smallest parameter value 

a, of the equation (py x ) x +%qy=0, and when — $■ — <n 2 , this in turn is (§ 1, 1) 

at least as great as the smallest of the equation y xa> +fa?y=0 under the same 
boundary conditions (16). Since for the latter equation ^=1 and Y 1 =c sin n%, 
it follows that % is greater than unity. Hence $ l t py x dx=% is greater than 
flqy 2 dos = l and D(y, 0, 1) is positive. 

In the non-definite case one must therefore have — : >7t 2 in at least a 

mm p 

part of the interval, and this assumption we shall now make. But we note 

here that in a sub-interval a 1} o^ such that y{a^) =y(a 2 ) =0, any value of q for 

which — F — <-, 7-5, makes the integral r„Hpyl — qy 2 )dx definite, as is 

mm P (Oji — «i) 2 Ja^rvx 
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readily shown by a comparison argument like that above. Hence, for given p 
and q, if throughout the interval 0, 1 the zeros can be taken close enough 
together (as e.g. may certainly be done where k has the same sign as k{%) and 
is taken large enough in absolute value, cf . § 1, I) , the integral for each sub- 
interval is definite, and hence D(y, 0,1) is positive. This is also true under 
the more general boundary conditions (2). 
On the other hand from the formula 

(pyl—qy 2 )dx=k Jcy 2 dx, (17) 

which is obtained by multiplying (15) by y and integrating under the condi- 
tions y(a 1 )=y(a 2 ) =0, we see that if k is positive and h negative in the 
interval, then D(y, a x , a 2 ) is negative; it is, however, still possible that for 
the larger interval 0, 1 the integral D(y, 0, 1) be positive. Were this the case 
(as will later be established provided k is large enough) we could deduce from 
Theorem IV and the formula 

D(y, 0, 1) Bs f\py 2 -qy 2 )dx=k Cktfdx, y(0) =y(l) = 0, (18) 

certain facts in regard to the movement of the subsequent zeros of particular 
solutions Y of (15) which vanish at #=0. Before proceeding to a detailed 
discussion of this matter let us compare some possibilities for the orthogonal, 
polar and non-definite cases. 

In the orthogonal case let us consider all k greater than some finite num- 
ber chosen less than the smallest characteristic (which may be positive or 
negative). As k increases, the subsequent zeros of any solution Y of (15) 
[7(0) =0] move to the left, new ones being added to the interval 0, 1. For 
the smallest value of k there are no zeros within the interval and when k 
passes through a parameter value k n the n-th zero enters. In the polar case 
there are two ranges of values of k ; for the range which extends from zero to 
+ co, a result precisely like that of the orthogonal case may be stated; for the 
range which extends from — co to zero, a decrease of k causes subsequent zeros 
to move to the left, there being no zero of Y(x, 0) present in the interval while 
the ti-th enters for k=k_ a . In the non-definite case we know that the results 
must be quite different and we may expect that the march of the zeros will not 
be monotone with A. In fact, there may be a range of values of k (L^k^L^) 
such that as X increases the number of zeros first decreases, then increases, 
then decreases and finally increases, the minimum number being a positive 
integer* 

* A reference to the proof of this last fact is given in the Introduction. 

38 
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As an example let us consider the equation 

y„+ [ (1007t) 2 + ^)2/=0, 2/(0) =i/(l) =0, 

where x is equal — 1 in the sub-interval 0, £ and equal 1 in the remainder of 
the interval. Such a function x(x) may be approximated by an analytic func- 
tion k(x) for which the corresponding equation has solutions with similar 
properties. "When Jl = the solution ?/ = sin IOOti* has 100 zeros in the interval. 
When %= (IOO71) 2 , we have the equation y xx =0 in 0, \ and ?/, ra . + 2(1007t) 2 ?/ = 
elsewhere, the solutions being respectively y = cx (c = const.), y = sin 100y27i#. 
There are are then approximately 50-\/2 zeros. For 2, = — (IOO71) 2 there are 
evidently the same number of zeros. More generally, when a, is in the interval 

— (IOO71) 2 , (IOO71) 2 the number of zeros is approximately 

£- [ V(1007t) 2 +^+ V(IOOti) 2 — X] , 

a function which has its maximum for X.=0 and decreases with increase of \"k\ . 
On the other hand when \%\ is greater than (1007t) 2 the number of zeros is 

approximately ^- V(100jt) 2 + \"k\, which increases with \%\. The minimum 

number of zeros is then approximately 50y2 and occurs for "k— (IOO71) 2 and 

— (1007t) 2 .* For values of n between 50 \/2 and 100 there are four parameters 
corresponding, to which exist solutions oscillating n times. 

Eeturning now to the discussion of the sign of D(y, 0, 1) we can prove 

Lemma I. If in an interval a, b, the function k is positive, except perhaps 
at the end points, then X can be taken so large that for solutions of (15), 

r b 
D(y,a,b)^\ (pyl—qy 2 )dx>0. 

"J a 

For, by taking A, large enough, q + Xk can be made as large as desired, 
except perhaps at the end points. Hence the zeros are as thickly strewn 
throughout the interval a, b as desired (§1, I). We note that in any sub- 
interval p\, /3 2 of a, b, in which there is a zero of y and for which the maximum 
of \y\ isJf, the minimum value of f^yldx is taken on when y is a straight line 

M 2 
joining (pi, 0) and (/3 2 , ±.M). This minimum value is -5 -=-. Since on the 

Pa— Pi 
other hand $ p e [qy 2 dx< (/? 2 — (3JM 2 max q, the value of D(y, (3 lt (3 2 ) is greater 

than M 2 \ — §- — ((3 2 — (3^ max q : hence when (/? 2 — p\) 2 is made less than 

Lp 2 — Pi J 

*By reference to Theorem IV it follows from this discussion that in the two intervals which are 
approximately — (1007r)*,0; (1007r) 2 , + 00 the value of J" 1 aiyHco must be positive or fiy*da><( l ,y 2 dx, while 
in the intervals which are approximately — 00, — (1007r) 2 ; 0, (100tt) s , f i y 2 d(v> C y 2 dw. 
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min p 

this integral is definite. It follows that since the first and last zeros of 

max q 

y(x) are as close to the ends of the interval a, b as is desired, those portions 

of the integral arising from these two end sub-intervals can be made positive. 

For any interval a 1} a 2 of oscillation [y{a x ) =«/(a 2 ) =0] it follows at once 

from (17) that when k and X are positive D(y, a lf a 2 ) >0. Combining these 

results we have the lemma. 

It follows in the same way that if k<0 and k ^£0 in an interval a, b f then 
'k can be taken so large negatively that for solutions of (15) the integral 
D(y, a, b) is positive. 

Coming now to an interval c, d in which k is negative, except perhaps at 
c and d, it is possible in an interval c-\-ri, d — q (?i>0 arbitrarily small) to 
take M' arbitrarily large and then choose a, so large that — (q + hk) >M' or 

\{py*)*\ = \-{<l+'Mi)y\>M'\y\, [{py x ) x >M'y if y>0]. (19) 

The arc y(x) may be shown to be sharply concave away from the *-axis.* It 
is possible to prove the following : 

Lemma II. If in an interval c, d the function k is negative, except per- 
haps at the end points, then /L can be taken so large that for a solution of (15) 

D(y, c, d)^ I (pyl—qy % )dx>0. 

This will be proved by showing that on taking 2, great enough it is 

V 2 
possible to insure that ^f is great at pleasure except perhaps in sub-intervals 

y 

whose length decreases indefinitely with — -. Since in the interval c + rj, d—v} 

the curve is concave upward for positive y and concave downward for negative 
y, there can not be more than one zero. The discussion may be separated 
into two parts according as y(x) has a zero or not. 

If y(y) =0, then without loss of generality it may be assumed that within 
the interval y, d—yj, y is positive and hence, since y m > 0, that y x is positive. 
Multiplying the inequality (19) by 2py x and integrating we get 

^VyAvyx)xdx>M' j 2ypy x dx, 
and hence 

p i {x)yl{x)>p i {x)yl{x)-p i (y)yl{ 7 )>M' [mm p} foyyjx 

— M'y 2 (x) mm p, (20) 

♦This is at once evident if the equation is taken in the form y x * -j- (q + X7c) y = to which it may be 
reduced by a transformation. 
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for all x in the interval y, d — v\. Similar reasoning establishes the same 

result in the interval c+v;, y. Hence, for the case that y vanishes, we have 

if M' min p 
throughout the interval c+>7, A— n the inequality Z\ > - 2 , which as 

y p 

noted above, is sufficient to establish the theorem. 

If on the other hand y has no zero in the interval, let us take the function 
positive and denote by y the point at which it takes on its minimum. It may 
be "assumed that y is not at the right-hand end of the interval; that special 
case could be treated in an analogous fashion. By multiplying y by a con- 
stant, y(y) can be made equal 1, while the sign of the integral D(y, c, d)' is 
not altered. Since y x (y)^0, we have y>l in the neighborhood of y; in fact, 
for any arbitrarily small but fixed r\ we can, by taking M' large enough, have 
y{y-\-Yi) >2. For, since y>l it follows from the inequality (19) that 

pyx>py*—pyAy)= I {py„).dx>M' \ ydx=M'(x—y). 

*J y %J y 

Integrating again, 

y x Ax> i py x dx>M' \ (x — y)Ax= ^— - — . 

y ms y %J y & 

MV 

Bv taking x— y> yi and choosing >1 we have y — 1>1 or y>2. From 

J o / == i o 2 max p 

this it follows that on multiplying the inequality (19) by 2py x and integrating 
we get by a process similar to that used in (20) the formula 

J-.x 2,M' 

2ypy x Ax>M'[y 2 (#)— 1] min p> —~y z (x) min p. 

if" 3M' min p 
Hence in the interval y+rj, d — r, we have— \ > -j % — . And since in the 

interval c+j?, y — v\ a similar inequality may be obtained, this completes the 
proof of the lemma. 

It follows in the same way that if k>0, except perhaps at the end points 
of the interval c, d, then % can be taken so large negatively that for solutions 
of (15) the integral D(y, c, d) is positive. 

We have then proved by these lemmas that in any sub-interval in which k 
has one sign, D(y) can be made positive by taking A, large enough positively 
and also by taking 2, large enough negatively. It follows that D(y, 0, 1) 
would be positive, and from the formula 

^ CkyHx = C\ py l-qf)dx^D(y, 0, 1) >0 

^0 «^0 
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and Theorem IV, that for | X | large enough the zeros of Y move to the left 
with increase of \%\. We can then enunciate 

Theobbm VII. There exists an integer n 2 such that for n>n 2 there are 
precisely two solutions of (15) (16) oscillating n times. 

If proper restrictions are imposed on G ($,%), a theorem similar to the 
preceding may he proved for the general equation 

(py x ) x +G(x,2 J )y=0. (10) 

The proof used in Theorems VI and VII, and Lemmas I and II is valid also 
under the following hypotheses on a function G'(x, /I) obtained by subtracting 
from G(x, A) a function q(x) • [G' (x, %)=G(x, 2,) — q(x)]. 

(a) In one or more sub-intervals of 0, 1, G' (x, 2,) is a monotone increasing 
function of X, such that lim G' (x, X)>0 and for at least a part of each sub- 

X=oo 

interval lim G(x, X) = + oo. 

X=co 

(/3) In the remaining sub-interval or sub-intervals of 0, 1, G'(x,X) is a 
decreasing function of % such that lim G'(x, /i)<0, and for at least a part of 

each sub-interval lim G' (x, X) = — oo. 

\=— CO 

(y) For values ?/, W of X there are n x zeros of the solution Y within the 
interval 0, 1 and for %>W, 

-£->4>iWO'(x,^), and for X<%",~ <-&(*)& (x, a), 

where <p 1} <£> 2 are positive functions denned in the intervals W, +co and — oo, A" 
respectively. 

Theoeem VIII. Under the conditions (a) (/3) (y) there exists an integer 
n 2 such that for n>n 2 there are precisely two solutions of (10) (16) oscillating 
n times. 

When the parameter value and corresponding solution of (15) (16) are 
complex, let us set hssa + n, y^u + iv. The differential equation resolves 
itself into the two following : 

(pu x ) x +qu-\-(jku — ikv=0, u(0)—u(l)=0, (21) 

(pv.).+qv+akv+i;ku = 0, «(0)=«(1)=0. (22) 

Corresponding to 2,^=0- — i% there is then a solution y^u—iv. On multiplying 

(21) and (22) in the first place by u and v, and in the second by v and — u 

respectively, adding and integrating, one obtains the two formulae 

f [p(ul+vl)— q(u 2 +v 2 )]dx = 0, f k(u*+v 2 )dx = 0. (23) 
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More generally one obtains by this process the formula 

[p(u x v—v x u)] x x zP — 'v J k(u 2 +v 2 )dx. (24) 

The orthogonality of two solutions (f]kY n Y m dx=0, m^n) holds as well 
for complex as for real characteristic numbers \, X m . Hence by separation 
into real and imaginary parts we get the 

Theoeem IX. If Y m (x, ^ m )=M m +w m , Y n =u n +iv n are two solutions of 
(15) (16), then 

C k(u m u n —v m v n )dx = 0, C k(u m v n +u n v m )dx = 0. 

Theorems of other types may be derived of which the following is an 
example : 

Theoeem X. 7/ in the interval 0, 1, the function k(x) changes sign once 
only, then the roots of the real and imaginary parts u n , v H of the solution 
Y n =u n + iv„ of (15) (16) separate one another. 

For, formula (24) may be written 



d' u 



p{u x v—v x u) = fk(u 2 + v 2 )dx, or p — ^-L — — C\(u 2 +v 2 ) 

J a dx V 2 Jo 



dx. 



From the hypothesis that k changes sign but once it follows that the integral 



d 



can not vanish within the interval, and since — - — - >0 except at the end 

dx 

at 

points, — is a monotone function and the theorem follows at once. We may 

note further that under the hypotheses of the theorem neither u nor v can 
gain or lose a zero at a point within the interval. For, were there a zero lost 
or gained at x=a, both y(a) and y x (a) would vanish, and formula (24) could 
be written f o k(u 2 + v 2 )dx=0, which would give a contradiction. 

§ 5. Reduction of the General Boundary Conditions to Normal Forms. 

In the preceding sections the discussion has dealt mainly with the simple 

boundary conditions y(0)y x (0)=y(l)y x (l)—0. To facilitate the discussion 

of the most general boundary conditions it is desirable to obtain normal 

forms.* The linear self-adjoint equation of the second order will be taken in 

the form (7t(x)u x (x)) x +T(x,X)u(x)=0, n(x)>0, (25) 

* The classification used in this section follows that of Hilbert and Haupt, loo. cit. The geometrical 
form into which the transformation is thrown was suggested by my colleague, Prof. H. P. Manning. 
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and if the linearly independent boundary conditions 

ajtt (0)+a 2 u 3l (0) + a s u (1) + a#» (1) =0, 
&«(0) +ftt*,(0) +ftw(l) +&u,(l) =0, 

are to be self-adjoint, the condition 



*(0) 



M(0) ^(0) 

5(0) 5,(0) 



■Ml) 



u(l) «,(1) 

i(i) 5,(i) 



(26) 



(27) 



must be imposed, where u, u are any functions with continuous derivatives 
which satisfy (26). We have the relation 

u(0) «,(0) _ a 3 tt 4 u(l) «,(1) 

5(0) 5,(0) _ & & 5(1) 5,(1) ' ( " ' 

as may be seen by applying the usual rule for multiplying determinants and 
noting from (26) that each of the four elements of one product determinant 
is the negative of the corresponding element of the other. Let us denote by 
B tj the determinant formed by taking the i-th and ;-th columns of the matrix 

a x a 2 a 3 a 4 

ft &. ft ft 

Considering the two determinants of the u's as the variables, we have from 
the theory of linear equations that a necessary and sufficient condition for the 
solution of (27) and (28) is that 

7t(l)B 12 -7t(0)B 3i =0. (29) 

Hence B 12 and B u are simultaneously zero or different from zero. 
Let us subject the dependent variable to the transformation 

u{x)=vi{x)y{x), yi(x)=f=0, (30) 

by which the number of zeros of the solution remains unaltered. The equation 
(25) takes on the self-adjoint form 

{Py*)z + G{%, ^)# = 0, p=7l>7 2 >0, G(X, ^,) =71)7^-1- 7t,57^ + r>7 2 , (31) 

and the conditions (26) are replaced by a similar set 



ri#(0) +y,y,(0) +y,y(l) +y 4 y,(l) =0, 
«iy(0) +5^,(0) +^2/(1) +8^,(1) =0. 

We have at once from (31) the important formula 



r x dT 2 , C x dG 2 , 
I ^- irdx— I ——y 2 dx. 



(32) 



(33) 
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Writing u lf u 2 , u 3 , u t , for u(0), u x (0), u(l), u x (l), and using a similar 
notation for the vj's and y's we can interpret our problem in relation to the 
tetraedron of reference as the investigation of the line (26) subject to the 
condition (27) or (29). If B 3i =0, the u 3 and w 4 can be eliminated from (26) 
at the same time, giving 

B w u 1 +B 2 aU 2 =0 or B u u 1 —B i2 ti 2 = 0, (34) 

either of these equations representing the plane determined by the given line 
and the edge 12(« 1 = 0, «. 2 =0). Hence the line (26) intersects this edge of the 
tetraedron. But if i? 34 =0 we have from (29) B 12 =0, and the straight line 
intersects also the edge 34. The equation of the plane determined by the 
given line and the edge 34 can be written 

B 13 u 3 +B u u 2 =0 or B 23 u 3 —B i2 u 2 = 0. (35) 

From (34) and (35) we see that under the hypothesis B 12 —B 3i = the boun- 
dary conditions may be written 

Case I. a , w 1 +M 2 =0, tm 3 +m 4 =0, (36) 

where 

•°1S n \i _ -"IS °i& 



#23 B i2 ' B u B 



42 



The parameters a and t may have zero or infinite values when the line (26) 
lies in a face of the tetraedron. 

In general, when any two B's with complementary indices are zero, the 
line intersects two opposite edges of the tetraedron and the conditions (26) 
may be reduced to a normal form similar to (36). 

The transformation (30) takes the form 

«i=>7i3/i> «2 = >722/i+>h2/2, u 3 =ri 3 y 3 , u i =Yi i y z -\-vny i , ^ =f= 0, *? 2 =£ 0, (37) 
and corresponds to a rotation of the face 2 of the tetraedron about the edge 
12, and of the face 4 about the edge 34, leaving the faces 1 and 2 and the 
edges 12, 13 and 34 unchanged. If we denote by A i} the determinants of the 
matrix of the coefficients y if B t of (32) we can read their values from the 
identity 



7\ 7t 73 7* 

§1 K $3 &4 



thus A 12 =Y[\B 12 , A 3i =v;lB 3i , etc. As may be seen from these formulae and 
(31), we have corresponding to (27) or (29), the relation 

p{l)A 12 -p(0)A 3i =0 (38) 

which, when added to (32), makes these boundary conditions self-adjoint. 
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The discussion of the boundary conditions under the transformation may 
be sub-divided as follows : 

I. If the given line intersects the edges 12 and 34 (B 8i =B n =0), the 
same will be true after the transformation since A u = J. 12 =0. To reduce the 
condition to the normal form (36) (which we shall in general use) the trans- 
formation is superfluous. This case is called the sturmian. 

From the formulae (37) it follows at once that the conditions (36) can by 
proper choice of j?i , j? 2 > % > Vi be reduced to one of the simpler forms 

(a) «/i = «/ 3 =0; (b) y 1 =y i =0; (c) y 2 =y 3 =0; (d) y 2 =y i =0, (39) 

these corresponding geometrically to the cases where the given line coincides 
respectively with the edge 13, the new edge 14, the new edge 23 or the new 
edge 24. 

II. If the line does not intersect the edges 12, 34, but does intersect 13, 
the transformation may be so determined that it will also intersect 24 ; e. g., 
by making the plane 2 pass through the intersection of the given line with the 
plane 4. In this case since m 4 is unaltered we have *? 3 =1, »74=0, and hence to 

equate -4i3— %(»7i#i3+>72-B23) to zero we need only to choose — = — =^-. After 

>7l -023 

the transformation we have A lS =A 2i —0, and by elimination in (32) we get by 
aid of (38) 

Case II. 2/i=%3, ^P(0)2/a=P(l)^, 



where 



h-- 



A 2S A 3i p(l)A 12 _ p(l)A 



23 



^12 4- u ' p(0)A u p(0)A Si ' 



III. If the line does not intersect any of the three edges 12, 34, 13, 
(B u ^0, #12 =£0, B 2i =f=0) we can determine the transformation so that the 
face 2 of the tetraedron shall pass through the intersection of the line and the 
face 3, while the face 4 passes through the intersection of the line and the 
face 1. Thus, in the new tetraedron, the line intersects the edges 23 and 14, 
and we have the relations A u =Yi z (ri 1 B u -{-YiJ3 ii ) —0, ^ 2 3=»7i('73#2s+>74-524) =0, 

yj yj 

which determine the ratios — , — desired in this case. The equations of the 
line can now be written 

Case III. y\=lp{^)yi, ip{^)yi——yz, 

where 



Z=- 



"■34 



A 4 



p(l)A 12 p{l)A 13 ~ p(0)As 2>(0)^34 

39 
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Theobem XL By a change of dependent variable the self-adjoint equation 
(25) with self-adjoint boundary conditions (26) (27) may be reduced to an 
equation (31) which is again self-adjoint, and for which the self-adjoint boun- 
dary condition may be written in one of three forms characterised respectively 
by Aw=A u =0; A 2S =A u =0; A u =A i2 =0: 

Case I. ay(0)-\-y x (0)=0, xy(l) -\-y x (l)=0, a, x constants (including oo), 
Case II. y(0)=hy(l), hp (Q)y x (0)=p (I) y x (l), h= constant, 

Gaselll. y(0)=lp(l)y x (l), Ip (0)y x (0) — — #(1), l=constant, 

The number of zeros of the solutions of the equation remains unaltered under 
this transformation. 

The value of the integral J -^— u 2 dx is the same as that of the corre- 
sponding integral ( -^— y 2 dx. 

Coeollaey I. In Case I the transformation may be so chosen that the 
boundary conditions can be written in one of the special forms 

(a) y(0) =y(l) = 0; (b) y(0) =y x (l)=0; 
(«) yAO)=y(l)=0; (d) y x (0)=y x (l)=0. 

§ 6. Oscillation Theorems for the Sturmian Boundary Conditions. 

We propose now to study oscillation theorems for solutions of the differ- 
ential equation under the first or sturmian case of the general boundary 
conditions (26) (27). As was shown in the preceding section this case, char- 
acterized by the relations B J2 =B u =0 (A 12 =A Si —0), may be reduced to a 
study of the equation 

(py,),+G(*,Vy=o, (41) 

under the boundary conditions 

'y(0)+O,(0)=0, Ty(l)+y„(l)=0, (42) 

where <r and t are constants, the important special forms 

(a) y(0)=y(l)=0; (b) y (0) =y x (l) =0; 1 

(c) y x (0)=y(l)=0; (d) y x (0)=y x (l)=0; J (4d) 

being considered as included for the values and oo. It was further shown 
that the boundary conditions may always be reduced to one of these special 
forms. 

We proceed now to a study of the nature of the dependence of /I on <t and t. 
If y and y are defined as the following fundamental solutions of (41), 

y<0,a,)=0,&(0,a)=l; y(0, X) =!,&(<>, A) =0, (44) 
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any solution y(x) can be written y = c 1 y-\-c 2 . For determination of c x and c 2 
we have on substitution in (42), 

<TC 2 + Cl = 0, t[c x £(1, a) + cj(l, X) ] +c^.(l, X) +<w.(l, a) = 0. 

A necessary and sufficient condition that there be values of c 1} c 2 satisfying the 

equations is that 

1 cr 

= 0. 



1 cr 



Hence each of the parameters cr, t is a one-valued function of the other and 
of Jl. If we assume that a zero of the solution moves continuously to the 
right or left, and note what happens as two consecutive zeros pass through 
the point x = 0, we see geometrically that for any given cr the condition 
ay(0) +2/a,(0) =0 must be satisfied at some stage of the process. It is possible 
by analysis to ascertain precisely what happens when cr or t passes over the 
range — oo, + oo. Corresponding to cr = — oo, 0, + oo one has respectively 
?/(0) =0, y x (0) =0, y(0) =0, and we shall show that under proper conditions 
one and only one zero of y(x) is lost or gained by the process. A similar 
change in % has a corresponding result. 

If in formula (13) we substitute from (42) and the formulae 

, 8y(0) da 9y.(0) _ dy(l) ^ dy^jl) 



obtained from (42) by differentiation with regard to \ we obtain the relation 



s. 



1 ™y*dx = -p(0)~y*(0) + p(l)^y'(l). (45) 



di 
Either of the parameters cr, t may be held fixed. If t is fixed — =0, and in 

order that a and cr be one-valued functions of one another over the range — oo, 

r 1 BG 
+ oo of cr it is then necessary and sufficient that the integral I ——y*d x have 

J OA 

one sign. From such considerations one may deduce a variety of theorems of 
which the following are examples. 

Theobem XIII. If a' m , a' m+1 denote two successive characteristic numbers 

of a solution of (41) for the boundary conditions 43(a), then in order that 

there be for all a an intermediate value of X corresponding to which there is 

precisely one solution of (41) for the boundary conditions ay(0) +^(0) =0, 

r 1 dG 
y(l) =0, it is necessary and sufficient that the integral J -x— y % dx have one 

sign for all solutions y{x) concerned. 
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Theorem XIV. // % m , X m+1 denote two successive characteristic numbers 
of a solution of (41) for the boundary conditions 

ay(0) +y x (0) = 0, 2/(1) =0, a — constant, 

then in order that there be for all t an intermediate value of /L corresponding 

to which there is precisely one solution of (41) (42), it is necessary and 

C 1 dG 
sufficient that the integral I -^-y^dx have one sign for all solutions y(x) 

concerned. 

Theorem XV. If corresponding to W,%" there are solutions of (41) (43a) 

with n x and n 2 zeros respectively, and if ( -^ — y 2 dx has one sign, then in the 

interval X, X' there is one and only one value of /I for which (41) (42) has a 
solution y{x) ivith n zeros (n 1 <n<n 2 ). 

We have seen in § 5 that by a transformation of the dependent variable 
the equation (41) may be thrown into the form 

— _ /*"' dG _ r a * dG 

(py x ) x +Gy=o, j a -fay idx =j a -fay 2(lx > p( x )>°> 

where the boundary conditions (42) assume one of the special forms 

y(a 1 )=y(a 2 )=0; y(a 1 )=y !l .(a 2 )=0; y x {a x ) =y(a 2 ) = ; ^(a^ =^(a 2 ) =0. 

As an extension of Theorems IV and IV A we have then the following: 

Theorem XVI. If a zero of one of the functions ey + y x , ty-\-y x is held 
fixed, then with increasing X the zeros of the other move closer to the fixed 

zero or further away- according as I -= — y 2 dx is positive or negative, the inte- 
gration extending over the interval between the zeros. 

From the preceding theorems we can deduce still further results. 

J" 1 dG 
-— — y 2 dx has one sign, then between any hvo con- 
o 3/1 

secutive Xs corresponding to solutions of (41) (43a) or of (41) (43d) there 

will be one each corresponding to solutions of (41) (43fo) and (41) (43c), and 

between tivo consecutive Xs for (41) (43&) or (41) (43c) there will be one 

each for (41) (43a) and (41) (43d). 

Corollary: Denoting by /L a and % d two adjacent characteristic values for 
solutions of the problem (41) (43a), (41) (43d) respectively, and by \,X C 
and X b , X c the next greater and next smaller sets for solutions of the problem 
(41) (43&), (41) (43c) respectively, then if the number of zeros of the solution 
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corresponding to % a is m, the number corresponding to X d is m — 1. When the 
integral is positive the number of zeros corresponding to \ and "k c is m, and 

to W h and X' c is m — 1, while if it is negative, the reverse is true. 

7)C 
If we consider the special form * of (41) where -^— > and Hypothesis A 

of § 2 is satisfied while G(x, — qo) <0 and = — oo in at least some portion of 
the interval, we can trace the various values of /I corresponding to successive 
suites of <r and t, and state 

Theobem XVIII. There exists under these hypotheses an infinite set of 
characteristic numbers A x < Jl 2 < /L 3 . . . ., with limit point at +qo only, for 
which exist solutions of (41) (42), the solution y n (x) corresponding to 
\(n — l, 2, 3, . . . .) having n — 1 zeros within the interval. 

The results for the orthogonal case (■§ 1, II) of the special equation 

(py,).+ (2+M0y=o, (46) 

are contained in the foregoing theorem. In the other cases it is preferable to 
investigate the zeros through a discussion of those of the equation obtained by 
the transformation used in proving Theorem XVI. We have 

(PV.).+ (q+U)y=0, p(x)>0, y(0)y x (0) =y(l)y x (l) = 0, 
(pyl — qy 2 )dx = | ky 2 dx = I ky 2 dx. 

If the new equation is of the polar form we see at once that there are 
precisely two solutions of (46) (42) with n zeros in the interval. If the equa- 
tion is of the non-definite form there are two integers n x , n 2 , such that for 
n<n x there are no solutions of (46) (42), for n>n x there are at least two, 
while for n>n 2 there are precisely two (Theorems VII, VII A). 

§ 7. Oscillation Theorems for Case II of the Boundary Conditions. 

We shall next consider the exceptional case B i2j ^a 4 ^ 2 — a 2 /? 4 = for which, 
as has been seen in § 5, the boundary conditions may be written 

y(0)=hy(l), hp(0)y x (0)=p(l)y x (l), ft=constant =£0. (47) 

Defining the particular solutions y(x, /I), y (x, "k) by (44), and substituting the 
solution y=c 1 y-\-c 2 y in (47), we obtain equations for c x and c % , 

c 2 =h[ Cl y{l,%)-\-c 2 y{l,k)}, ^(0)c 1 =p(l)[c 1 ^(l,^)+c 2 f,(l,X)]. (48) 
A necessary condition for a solution is the equation 

hy(l,X) hy{l,X)-l 

p(l)y x (l,k)-hp(0) p(l)y x (l,V 

♦This is a slightly more general condition than that imposed by Birkhoff. Loo. cit. 



D= 



= 0. (49) 
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Since for the two solutions we have the well-known formula 

P(%) [%(%> %)y(x, ^)—yA x > *>)V(x> *) ] = constant = — p(0), (50) 
this condition reduces to 

h'p(0)y (1, Jl) -2hp(0) +p(l)y x (l, I) =0, 

from which we get by again using (50) 

i ± VF(i)V-y,(i>^)y(i,^) (51) 

y(l,^) Vp(0)y(l,X) 

So long as ^(l, ^) and «/(l, X) have opposite signs there will be two values of 
h for each value of % ; when they have the same sign there will be none. We 
may distinguish three critical cases for which the solutions pass from real to 
complex : 

(a) If y(l, /I) =0, y x (l, A,)=£0, the equation (49) can by means of (50) 
be written [hy(l,X)—l] 2 p(l)y x (l,X)—0, and since y and y x can not vanish 
together, it follows that hy(l, /L) — 1 = 0. Hence every element of the deter- 
minant D vanishes except that in the lower right-hand corner; it follows then 
that c 2 = 0. The solution y{x)z=c-${l,X) vanishes at x=0, and since h=f=0, 
we have from (47) that it vanishes also at x=l. The function y(x) is there- 
fore a solution of (41) (43a). 

(/3) If y(l,/i)^=0, y x {l,X)=0, we have in the same way ^=0 and the 
solution y=c 2 y also a solution of (41) (43d). 

(y) If y(l, %) =0, y x (l, /I) =0 it follows immediately by the same reason- 
ing that all elements of D vanish. Hence c x and c 2 may take on any values 
and /I is a double characteristic number. 

On multiplying equation (41) by y and integrating under the boundary 
conditions (47) we obtain the relation 

f pyldx = C Gy 2 dx. (52) 

Jo J 

By substitution in (13) from formulae (47) and the further formulae, 

^ = ^-1^. i»<0)*<0> + *<0)*#U(l>^. («> 

obtained from (47) by differentiation with regard to %, we obtain the funda- 
mental formula 

r f 1 3G 2J , 2 r 1 dG 2 , r l dG 2,7 

dh U-aa^ h J w ydx I d% ydx 4) 

tU 2p(0)y(0)y x (0) 2p{l)y{0)y x {l) "2p(0)y(l)y,(0) ' 
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Let us denote by \ a characteristic number of our problem for h = 0; this 
is a solution of the sturmian case (41) (43b). Near it, as we know from 
Theorem XVII and its corollary, there is another characteristic number a c 
corresponding to h — as and to (41) (43c), the two solutions for A b and a c 
having the same number of zeros. The parameter value A b may be equal to, 
greater than, or less than a c . Of the aggregate of characteristic numbers for 
the two problems (41) (43a), (41) (43d) let us denote by A M the greatest of 
those smaller than ?. 6 (or a c ) and by X ad the smallest of those greater. Within 
the interval A ad , X ad it follows from Theorem XVII that there is but one zero 
and one infinity of h ; in other words, but one solution for each of the prob- 

-— — y 2 dx has one sign 
o d/L 

we are now in a position to prove that in this interval of a the function h(X) 

defined by (51) is monotone on both of its branches (the one including h = 

and the other h=<x>). 

To prove this let us in the first place note that for h = or h — cc we have 

a sturmian case for which existence theorems have already been established, 

dh 
and then show that as h passes through either of these values, -=- does not 

aA 

change sign. Wesee at once from (47) that y x (l,h{0)) = 0; hence ?/(l,/l(0))=jk0, 

and it follows from (47) that y(0, A{h) ) changes sign with h. A reference to 

dh 
the first part of (54) shows that -^ will then retain its sign. In the same 

a A 

way y(l, M 00 )) =0, y x (l, X(oo))=fc0; y x (0, A(h)) changes sign as h goes 

dh 
through infinity and -j- retains its sign. This establishes the result since it 

dA 

follows from (54) that while j -- — y 2 dx retains one sign the only possibility 

Jo OA 

dh 
of -rr changing sign is when h or y(0) or y x (0) changes sign. 
a a 

By definition y(0, Jl(0)) =0, y x (l, a(0)) =0, and without loss of generality 

it may be assumed that y x (0,A(0)) >0. The discussion then divides itself 

into two parts according to 

Hypothesis I. y(l, a(0) ) >0. Hypothesis II. y(l,A(0))<0. 

Under the first hypothesis the number of zeros of y (including that at 

a;=0) is even, while under the other it is odd. Roughly speaking, we shall see 

that a gain or loss of a zero comes when h goes through or go. Concerning 

the function A(h) there is now sufficient data to sketch the graph. We shall 

i d/"y 
first discuss the problem for the assumption j -^ — y 2 dx>0 and Hypothesis 1. 
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dh 
It has been proved above that j- can not change sign on either branch of 

dX 

the function. On the branch containing h=0, A=A b we see from (54) that 

dh 

-rr <0, since for that particular point y(l) >0, ^(0) >0. For h = <x>, a—a c , 

aA 

the number of zeros is the same as that for h=0 (Corollary, Theorem XVII) 

and hence under Hypothesis I, y x (l) and y(0) have opposite signs. From the 

dh 
second part of (54) we see then that - 7 ->0on this other branch. Since at the 

dA 

ends of the interval h is not 0, we know that* in cases (a) and (/?) either y(0) 

or y x (0) will vanish according as A ad , X ad belong to the problem (41) (43a) or 

dh 
(41) (43d) ; it follows then from (54) that ^ =qo. 

u/A 

Under Hypothesis II it is readily shown by the same processes that in an 

interval which we shall call \ d , ?J ad to distinguish it from the other, the branch 

of the function A(h), which contains h = 0, is monotone increasing and the branch 

containing h= oo is monotone decreasing! As a increases this form of curve will 

i d/~* 
always alternate with that obtained under Hypothesis I. When I ~^~ y 2 dx<0, 

Jf) OA 

the second form of curve occurs under Hypothesis I, and the first under 
Hypothesis II. 

Since by definition the intervals can not overlap, the curves can not 
overlap. If the upper bound of one interval is the lower bound of the next, 
we have a double point for a [cf. (y) above], and one branch of one curve 
unites with one of the other to form a function monotone throughout. This is 
what takes place, for example, in the case of the solutions of the equation t 

*In case (7), 3/(0) and j/»(0) may be chosen arbitrarily; cf. next succeeding foot-note. 
f For this special equation y = ry= sin V\a>, y ^ cos VXaj, and equation (49) becomes 

2h — (W+l) cos Vx = 0, 

from which we obtain the formula ft = sec Vx±tan Vx. For every positive X the function h is double- 
valued. The interval X« s , X'«s is (2m)" it-, (2m + l) 2 7r 2 , while at (2m + 1/2) V, h becomes zero on one branch 
and infinite on the other. The monotone decreasing branch 7j x = sec VX — tan VX may be considered as 
joined at both ends to monotone decreasing branches in the next intervals. These intervals (2m — 1) 2 jt, 
(2»t) 2 7r 2 and (2m + l) 2 7r 2 , (2m + 2) 2 jr 2 are of the type X afl ,X'„,j. The function h t decreases from +co to — 00 
in the interval (2m — 1/2) V, (2m + 3/2) V. The function h 2 = sec Vx + tan Vx increases monotonely 
from — 00 to +00 in the interval (m-\-l/2)-ir 2 , (m + 5/2) s ir, the branches of adjacent intervals uniting 
as in the other case. Each of the curves \ and h. 2 cuts two of the other set orthogonally, the points of 
intersection occurring at the end points m 2 jr 2 of the interval X 0( ( , \',„, . 
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It is bow easy to write down oscillation theorems for this case. We note 
in the first place that under Hypothesis I the number of zeros is increased by 
unity as h goes through zero, and under Hypothesis II it is decreased by unity 
so that for h>0 the number of zeros is always even, and for h<0 it is odd. 
There are then two solutions with an even number of zeros in the one case and 
two with an odd number in the other. 

Theorem XIX. If in an interval of X there exist two integers m x , m 2 
positive or zero such that there are solutions of 

(py x ) x +G(x,2,)y = (55) 

and (43a) with m 1 and m. 2 zeros respectively within the interval 0, 1, then 

1 0/~f 

under the hypothesis that S ~^—y 2 dx has one sign throughout, there are for 

Jo CM, 

h>0 two solutions of (55) (47) when m is even, and none when m is odd 
(m 1 <m<m 2 ) : for h<0 there are two solutions when m is odd aud none when 
m is even. 

The special case where -~— >0 and G becomes negative for Jl= — oo, and 

moreover negatively infinite in at least a part of the interval and G becomes 
positively infinite for /l= + oo in at least a part of the interval, includes the 
case discussed by Birkhoff and the detailed theorems derived by him hold also 
here.* 

The orthogonal problem (§ 1, II) of the equation 

{py,),+ (Q+M)y=o (56) 

is contained in the special case just discussed. In discussing the polar case 
we note that it follows from the special case of the formula (52), 

0< I (pyl— qy l )dx—% C Jcy 2 dx, (57) 

Jo Jo 

r 1 dG r 1 

that j ——y 2 dx^ I ky 2 dx has the same sign as %,. Hence /l=0 is not included 

in the range of values (unless one includes the solution y=0). 

Theokem XX. In the polar case of equation (56) there are precisely two 
solutions satisfying the boundary conditions (47) and oscillating n times 

(» = 0, 1,2, ....)• 

In the non-definite case we have from Theorem VII that for sufficiently 
large values of |/l| the integral on the left of (57) is positive and hence we 
can state 

* Case II corresponds to III, p. 269 in Birkhoff's article, loo. cit. 
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Theorem XXL For the non-definite case of equation (56) there exist 
two integers n x , n 2 (%>%^0) such that for n<n x there is no solution of 
(56) (47) with n zeros in the interval; for n>n r there are at least two, while 
for n>n 2 there are precisely two. 

§ 8. Oscillation Theorems for Case HI of the Boundary Conditions. 

There remains the most general of the three normal forms obtained for 
the boundary conditions, viz. : 

y(0)=lp(l)y x (l), lp{0)y.(0)=-y(l), Z = constant =£ 0. (58) 

Since the discussion follows the same lines as that of the preceding section it 
will be abbreviated. If we define the two particular solutions y(x, a), y(x, A.) 
as in (44), we get in place of formulae (49), (51), 

lp(l)y x (l,X) lp(l)y,(l,*,)-l 

y(l,l)+lp(Q) y(l,l) 

(59) 



D' = 



= 0, 



t _ 1 V&(1, *,)£(!, Jl) 



p{l)y x {l,X) Vp(0)p(l)y,(l,Jl) 

The critical values are when y x (l, 1) =0 and y(l, /t) =0 and the cases may be 
classified as before. 

(a) If 1/(1, 2,) =0, y x (l,7i)^=0, equation (59) can be written 

j7(l,a)[Zip(l)^(l,a)-l] 2 = or lp(l)y x (l,l)=l. 

It follows that c 1 = and the solution y(x, %)=c 2 y(x, "k) is a solution of the 
equation with the sturmian boundary condition (43c). 

((3) If y(l, /t)=fc0, y x (l,X)=0 we have c 2 — and y^=c-$, which is a 
solution of (41) (43b). 

(j/) If 1/(1, a)=0, y x (l, h)=0, it follows as before that D' vanishes 
identically and ?„ is a double parameter value. 

To replace formula (53) and (54) we have 

di 0. gx ydx l J. w^ = "J. a* ^g 

<U 2p(0)y.(0)y(0) 2y(l)y(0) ~ 2p(0)y.(0)p(l)2/.(l) " ( ' 
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ji 
It is readily shown as in § 7 that -^- does not change sign as I goes through 

the values or oo. Let us denote by \, A d two adjacent characteristic num- 
bers for the cases 1 = and l = co, respectively (in other words for (41) (43a), 
(41) (43d) ), by \ e the greatest of the aggregate of characteristic numbers for 
(41) (43b) and (41) (43c) which are smaller than \, and by X ie the smallest 
of the aggregate larger. It is readily shown as in the previous section that 
within the interval \„ , X' be , I is a monotone function on the branch containing 
/ = and on the branch containing h=co. Let us consider first the case 1=0; 
then i/(0, A(0))=0, y(l,A(0))=0 and we can assume y x (0, A, (0) ) >0. There 
will be two cases to distinguish according as we make 

Hypothesis I. y x (l, a(0) ) >0 ; Hypothesis II. y x {l,A\{0) ) <0. 

Under Hypothesis I the number of zeros is always odd, and when I goes 
through zero from negative to positive it may be seen from (58) that two 
zeros of the solution y{oc) are lost. Under Hypothesis II the number of zeros 
is even, and, as I increases through zero, two zeros of y(x) are gained. It fol- 
lows from the last part of (60) that the integral I -~r~ y 2 dx and — ~- - have 

opposite signs. This fixes the sign of -=^ on one branch. To fix the sign of 

-77- on the other branch let us consider ^(oo)^^. Under Hypothesis I we can 

Q/At 

argue from the corollary to Theorem XVII that y(0, A(°°)) and t/(l, A(oo)) have 

the same sign, and from the third part of (60) that when | -^— yHx > 0, — )^1L 

Jo OAj dX 

is positive. Under Hypothesis II we can easily prove that the situation is 
reversed, -r- being positive on the branch through h=0, and negative on the 

other. "When I -5^- y 2 doc < the two types of curves are interchanged. With 

•'O OK 

increasing a we see then that in both cases curves of much the same form as 
the two varieties in the preceding section alternate with one another. If 
adjacent intervals have the same end-point, a branch from the one will unite 
with a branch of the other.* 

From these data various theorems may be deduced, of which the following 
is typical. 

♦For the special equation y x ,-\-\yz=0 one can set up the formula for l(\) as in foot-note, p. 312. 
It is found that I = — esc V\ ± cot Vx/VX and that these monotone functions extend from — <x to -f-oo , 
each of the one set cutting two of the other set in double points X= (n-f-l/2)V of l(\). 
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Theobbm XXII. If in an interval of 1 there exist two integers faK^, 
positive or zero, such that there are solutions of (55) (43a) with fa and fa 
zeros respectively within the interval 0, 1, then 

r 1 dG 

(1) Under the hypothesis that \ -~r-y 2 dx be positive, there are, when I 

is positive, two solutions of (55) (58) which have 2m or 2w — 1 zeros in the 
interval ( l u 1 <2m — l<2m<fa) and there are, when I is negative, two solutions 
which have 2m or 2m-\-l zeros (fa<2m <2m-\-l < fa) . 

1 CI/"* 

(2) Under the hypothesis that j -^r-y 2 dx be negative, there are, when I 

is positive, two solutions of (55) (58) which have 2m or 2w + l zeros in the 
interval (fa<2m<2m-\-l<fa) and there are, when I is negative, ttvo solutions 
which have 2m or 2m — 1 zeros (fa<2m — l<2m<fa). 

One can make the results of this theorem more specific by giving the con- 
ditions necessary to characterize the branch of the function Z(/l) which is 
involved. Birkhoff has done this for the special case treated by him (loc. cit., 
p. 269, I, II), but we shall content ourselves with stating that the same classi- 
fication may be made in the general case treated here. 

Theorems for the special equation (56) analogous to Theorems XX, XXI 7 
can be at once written down. 

Brown University, January, 1918. 



